Local properties and Density of States in the two-dimensional p-d Model 
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The p-d model which well describes the CUO2 planes of the high-Tc superconductors is studied 
by means of the Composite Operator Method (COM). The relevant quasi-particle excitations are 
represented by composite operators. As a result of taking into account spin excitations we find a 
p— like band near the Fermi level. The dispersion of this band gives a Fermi surface which is in 
good agreement with the experimental measurements. Due to the strong mixing of the relevant 
excitations, the spectral weight of this band is reduced and gives a large Fermi surface in the 
moderately doped region. The dependence of the calculated physical quantities on model parameters, 
temperature and doping, is in a very good agreement with the available Quantum Monte Carlo 
results. 
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I. INTRODUCTION 

Many unexplained exotic features of transition metal 
oxides containing Cu02 planes have inspired the 

developing of various theoretical methods to solve, at 
least approximately, the models which describe such com- 
pounds. 

In the last years, a big effort has been devoted to ex- 
plain experimental results in terms of Hubbard-like or 
t-J models and many interesting results have been suc- 
cessfully achieved. The Hubbard model was suggested 
as low energy derivative of the Emery (or p-d) three 
band model |^ . The Emery model takes into account the 
crystalline structure of Cu02 planes and initially retains 
the orbitals of p and d electrons at oxygen and copper 
sites, the d orbital giving the strong Coulomb repulsion 
when occupied by two electrons. Due to the more de- 
tailed structure of the p-d model, one can expect that 
its solutions can give us a more rich physical picture of 
cuprates and can explain a wider range of experimental 
results than solutions of Hubbard-like or t-J models. 

It is believed that the strong correlation of electrons in 
the copper d orbitals leads both to the insulating antifer- 
romagnetic state at low density of doping carriers and to 
superconductivity at higher doping. The doping depen- 
dence of the electronic state is not rigid band-like. This 
suggests that the doped hole changes the state of nearby 
localized copper d electron and reduces the energy of the 
system if compared to the case of rigid band doping. 

A detailed mechanism of binding d electron to the 
doped hole was proposed by Zhang and Rice As 
a result of binding, a local singlet forms and it can move 
through the lattice. Further theoretical investigations of 
the band dispersion of such singlet in Cu02 planes gave 
a successful description of ARPES results 

However, to recover the doping and temperature de- 



pendence of micro- and macro-scopic observables, one 
needs to study not only such singlet excitations, although 
their energy lies near the Fermi level, but also the wider 
basis of excitations which are possible in the plane. 

Recently, a technique called the Composite Operator 
Method (COM) has been proposed |l2|-|l|] to describe 
the local and itinerant properties of strongly correlated 
systems. A physical model is usually described in terms 
of elementary particles and some interaction. However, 
at level of observation the identity of the original par- 
ticles is lost; the macroscopic behavior of the system is 
described in terms of new excitation modes. When the 
interaction is strong the properties of these excitations 
will be very different from those of the original particles 
and hardly obtainable by a perturbation scheme. By fol- 
lowing this scheme, in the COM a set of field operators 
is taken as the basis in which the theory is developed. 
These fields are chosen in order to describe the exper- 
imentally observed properties and are called composite 
fields since they are constructed from the initial set of 
particles. The properties of these composite fields are 
dynamically determined by the interaction and by the 
boundary conditions and must be self consistently calcu- 
lated. In this process a special attention is put to pre- 
serve the symmetry properties of the model [|l5|,|l6| . Par- 
ticularly, the conservation of the Pauli principles plays 
a fundamental role. As discussed in Ref. although 
the method has many points in common with the pro- 
jection method | [l7|] , the method of equation of motion 
llTsI and the spectral density approach jl^], remarkable 
differences arise according to the different schemes of self 
consistency. The method has been applied to the study 
of several models, like Hubbard t-J (2l) and 

p-d fllllll models. 

The insufficiency of studying only those electronic 
states with energy close to the Fermi level was shown 
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in Ref. [|I2[. The reduced p-d model, where the transi- 
tions to the lower Hubbard level are neglected, was stud- 
ied within a four-pole expansion for the Green's function. 
The fields were taken as the bare p electron, the upper 
Hubbard field, the p electron field with spin-flips at cop- 
per sites and the p electron operator with p-d charge- 
transfer. 

In this work we investigate the full p-d model with 
a four-component field constituted by p electrons, lower 
and upper Hubbard operators for d electrons and a fourth 
field which describes p electrons dressed by spin fluctua- 
tions of d electrons. Some preliminary results were given 
in Ref. p^ , where in particular we showed the neces- 
sity of introducing the fourth field. The model treated 
in terms of only the first three fields gives results not in 
agreement with the numerical simulation data and that 
do not give a satisfactory description of the experimental 
situation. 

In Sec. |ll| we present the p-d Hamiltonian and intro- 
duce the excitation modes represented in terms of com- 
posite field operators. The Green's function for these op- 
erators is calculated in the four-pole approximation and 
depends on a set of parameters expressed as equal-time 
correlation functions. The study of these parameters is 
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given in Sec. HI, where by using the symmetry imposed 
by the Pauli principle we derive a set of equations which 
allow us to obtain a complete self-consistent solution. 

The developing of numerical methods, as Quantum 
Monte Carlo (QMC) and exact diagonalization ||2^, of- 
fers a new challenge to theoreticians. Before looking at 
the physical implication of the model one should be confi- 
dent with the solution by calculating some quantities and 
comparing them with the data of numerical simulations. 
According to this, in Sec. p^V A we make a detailed com- 
parison of our results with the available data obtained 
by QMC studies. The comparison shows a very good 
agreement signaling us that the choice of the composite 
field and the self-consistent procedure adopted are able 
to give a reason able s oluti on of the model. We then pro- 
ceed in sections IV B and [V C to analyze some physical 
properties. Section V is devoted to concluding remarks. 
Details of calculations are given in Appendix. 



II. THE P-D MODEL WITHIN THE 
FRAMEWORK OF COMPOSITE OPERATOR 
METHOD 

In this section we consider the p-d model and we 
present the relevant features of the Composite Opera- 
tor Method. Starting from the tight-binding model ||] 
composed of d electrons in 3dx2-y2 copper orbitals and 
p electrons in 2px^y oxygen orbitals in a Cu02 plane and 
considering only bonding p electrons p^ , we keep the 
terms describing the strong intra-atomic Coulomb repul- 
sion at Cu sites and the p and d electrons hopping. 

Therefore, the original Hamiltonian of the p-d model 
can be written as 



E 



{Sd - fi)d'^ {i)d{i) + (£p - fi)p'' {i)p{i) 



(2.1) 



where in the spinor notation and d\{i) create a p 

electron and a d electron on site i = (Ri,ii) with spin 
projection a G {t, |}, respectively, /i is the chemical po- 
tential, the U term describes the Coulomb repulsion at 
Cu sites, ncr(i) = d\[i)d„{i) is the charge density opera- 
tor of d electrons with spin a. The t term describes the 
p-d hopping and p''{i) is defined by 

p'{i)^Y.^^,p{J)■ (2.2) 
3 

Here 7ij is given by 

7». = 7^ / d^ke^'^-(^-^^)7(k), (2.3) 



where fis is the volume of the first Brillouin zone, 7(k) = 
yjl — a(k) and where for a two-dimensional quadratic 
lattice with lattice constant a 



a(k) — i [cos(fc:Ea) -I- cos(fcya)] 



(2.4) 



Due to the strong correlation among electrons we in- 
troduce the composite field 



vE'(z) 




where 



^(z) = [l-n(z)]d(i), = n(i)d{i) 



(2.5) 



(2.6) 



are the Hubbard operators which describe the basic ex- 
citations n (i) = <-> n (i) = 1 and n(i) = 1 ^ n{i) — 2 
on the lattice site i, respectively. The fourth field is cho- 
sen as 



Ps{i) = (7knk{i)p^{i) 



3c 



l22 



3b 



^33 



(2.7) 



The parameters b and c and the quantities I22 and 133 
are given in Appendix. 

We have introduced the charge- (/i = 0) and spin- 
(/i = 1, 2, 3) density operator of d electrons 

nf,{i) = d^i)a^dii), (2.8) 

and we are using the following notation 

a^ = (l,?), a^ = {-l,<j) (2.9) 

with cTfe (fc = 1, 2, 3) being the Pauli matrices. 
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The choice of the composite field (2.5) is dictated by 
the following considerations. The strong intra-atomic 
Coulomb interaction at Cu sites induces a splitting of 
the d band into the lower and the upper Hubbard sub- 
bands. The large covalence between oxygen and copper 
electrons leads to a large fluctuation of the energy of p 
electrons at O sites. The field ps, describing p— elec- 
tronic excitations accompanied by the nearest neighbor 
d— electron spin fiuctuations, represents electronic exci- 
tations associated with the Cu-0 bonds. 

The Heisenberg equations of motion for the composite 
field ^' {i) are 



■9 w•^ 



(£d-A^)e(^) + 2^pT(^) + 2^7^(^) 
{e.p - fJ.)ps{i) + £ppP'^{i) + £p!i^{i 



(2.10) 



pTj^ ^p and the 



where the coupling constants £pp, Sp^, & 
higher order operators 7r(i) and Ks{i) are defined by the 
following relations: 



-pp 



etc 
I22 ' 



Spi 



3c(ep - Ed) 



'22 



(2.11) 



{{SJ{i),^\j)}E.T.)=0. 



(2.17) 



In the framework of the pole approximation we ne- 
glect the nonlinear term dJ{i) in the equation of motion. 
Then, for a translational invariant system the Fourier 
transform S'(k, uj) of the retarded propagator has the fol- 
lowing expression 



5(k,c.) = 



1 



UJ - e(k) 



(2.18) 



where £(k), the Fourier transform of the energy matrix, 
can be expressed as 



£(k) = m(k)/-i(k). 



(2.19) 



In order to simplify the notation we introduced the so- 
called normalization matrix /(k) and the so-called mass 
matrix r7i(k) 



I{k)^F.T.{{^{t),^Hj)}E.T.) 



m{k)=F.T.({i-^it),'f\j)}E.T. 



(2.20) 



(2.21) 



3b{ep - £,,) 



6t 



^33 



b 

I33 



c 
I22 



(2.12) 



(2.13) 



Ks{i) = akd\i)akP^ {i)p^ [i) 

-akp'^'\i)akd{i)p'^ {i) + aknk{i)d^ («)• (2-14) 
Let us introduce the thermal retarded Green's function 
S{i,j)^{R[^m\j)\), (2.15) 

where R is the usual retarded operator and the bracket 
(• ■ •) denotes the thermal average on t he gr and canonical 
ensemble. Use of equation of motion ( 2.10| ) will generate 
an infinite hierarchy of coupled equations and we need 
to introduce some approximation . To this purpose, let 
us split the Heisenberg equation (2.10) into a linear and 
nonlinear part 



d 



-^,{^)=Y^e{^,m{J) + 5J{^), 



(2.16) 



where the equal-time correlation matrix e{i,j), the so- 
called energy matrix, is fixed by the requirement that 
the nonli near term SJ{i) is orthogonal to the fundamen- 
tal basis dU): 



where the symbol F.T. denotes the Fourier transform. 
Straightforward calculations give the expressions of these 
matrices. The results for a paramagnetic state are given 
in Appendix. The Green's function ( 2.18| ) can be put in 
the spectral form 



^(k,^) = E 



r(») 



(k) 



(2.22) 



The energy-spectra £'„(k) are the eigenvalues of the 
£(k) matrix and the spectral functions cr^"^(k) are given 

by 



J") 

'ab 



(k) = ^A,„(k)A-i(k)/,6(k), 



(2.23) 



where a, b — 1, ... ,4 and where the columns of the A(k) 
matrix are the eigenvectors of the £(k) matrix. 

Summarizing, in our scheme of calculations firstly we 
choose a fundamental basis of Heisenberg composite op- 
erators. The propagator for this basic field is evaluated 
in a pole approximation, where the incoherent part is 
neglected. The main ingredient of the calculation is the 
energy matrix. The knowledge of this function will al- 
low us to calculate the excitations of the system and the 
spectral functions. The entire process requires a self- 
consistent procedure which will be discussed in the next 
Section. 
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III. CORRELATION FUNCTIONS AND 
SELF-CONSISTENT EQUATIONS 

In this section we introduce the correlation functions 
and give the expressions of the self-consistent equations 
needed to calculate the fermionic propagator. 

As shown in Appendix, the Green's function depends 
on a set of parameters expressed as static correlation 
functions of composite opera tors . Some of these oper- 
ators belong to the basic set (2^) and their expectation 
values are directly connected to matrix elements of the 
Green's function. Other operators are composite fields of 
higher order, out of the basis (2.5), and their correlation 
functions must be evaluated by other means. This aspect 
is a general property of the Green's function method [ p7[ . 
These functions refer to a specific choice of the Hilbert 
space and one must specify the proper representation. In 
particular, the states must be constructed in such a way 
that the relations imposed by the Pauli principle must 
be satisfied at level of expectation values. 

Let us introduce the correlation function 



(3.1) 



By means of the spectral theorem and by recalling Eq. 
( 2.22 ), this function has the expression 



X [l + r„(k)]a(")(k), (3.2) 

where 



T„(k) = tanh 



2kBT 



(3.3) 



The parameters directly connected to the Green's func- 
tion are: n^, a^, &, 6s, c, D. By means of the definitions 
given in Appendix, we have the self-consistent equations 



Ud = 2(1 - C22 - C33), 



3r 

"-'14 + 1~*-'12 
122 



6 = 



36 



^33 



•-^13' 



6s — C^4 



13' 



3r 

34 + ~r~ ^^22 
122 



(3.4) 
(3.5) 
(3.6) 



23J 



'~r~ [^23 + ^"^3] 

J. 33 



^12^ 



D = I: 



33 



33- 



(3.7) 

(3.8) 
(3.9) 



We are using the following notation 



(3.10) 
(3.11) 

(3.12) 



with 



«^. = 7^/ rf^ke^'^(^-^^)«(k). (3.13) 

The parameters not directly connected to the Green's 
function are d^, /, Xs, defined in Appendix. In princi- 
ple there are different ways to calculate these parameters: 
decoupling approximation, projection on the basis, use of 
linearized equations of motion. However, as discussed in 
Rcfs. , there is only one definite way to fix them: in 

order to have the right representation the Green's func- 
tion must satisfy the equation 



hm 5(z,j) = (*(*)*'(*)), 



(3.14) 



where the r.h.s. is calculated by means of the Pauli prin- 
ciple. 

Use of this equation leads to the following self- 
consistent equations 



23 



0, 



24 



_ 0,^7 0Oi2'-'22 
— t><^Y2 J • 

122 



C: 



34 



-^13'-'33 
-^33 



TIT = 4 - 2(Cii + C22 + C33), 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



where ut = np + nd is the total particle n umb e r an d 
i s the number of p electrons. Equations (3.4)-(3.£) and 
(|3.15| )- (|3.18|) constitute a set of coupled equations which 
will fix the parameters in a self-consistent way. 

Summarizing, we have ten parameters and ten self- 
consistency equations which allow us to compute the 
Green's function and the properties of the system in a 
fully self-consistent way. 

Results of calculations are presented in the next Sec- 
tion. As as a comparison we also present results in the 
case where we do not take into account the Pauli princi- 
ple and we express the parameters ds, /, Xs by means of 
the following decoupling equations 
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2C23 



Xs 



'2 [C22 



2C23 



^33! 



(3.19) 



(3.20) 



(3.21) 



It should be noted that for certain values of the exter- 
nal parameters some instabilities appear in the iteration 
procedure and no solution of the self-consistent equations 
is found. This aspect is not present when the decoupling 
scheme ( ^.19 )-(3.21) is used and may be related to the 
fact that more asymptotic fields are necessary in order 
to have a proper representation for the Green's function 



IV. RESULTS 

This Section is organized as follows. The first part 
is devoted to compare the results of our calculations 
with the data of numerical analysis by Quantum Monte 
Carl o and exact diagonalization. As it will be shown in 
Sec. IV A, for all local quantities the agreement between 
numerical and COM results is excellent. Once we are 
confident to have a reasonable solution of the model, we 
go to the next step where we study some physical prop- 
erties in order to verify if the p-d model is a realistic 
model for cuprate su percon duc tors. This is done in the 
second part, Sections IV B and IV C , where the density 
of states and Fermi surface are analyzed. In this study 
the values of the model parameters have been taken ac- 
cording to the results suggested by ab initio calculation 
|8|: i = leV, U = 6eV, A = 2eV, where A = Sd + U-Ep 
is the charge-transfer energy. 

In this work all energy are given in units of t and mea- 
sured with respect to the atomic level Sp = 0. 



A. Comparison with QMC data 

Here we compare COM results with QMC calculations 
and Lanczos diagonalization results j3^], | ]3l[ |. 

We introduce the squared local magnetic moment for 
d electrons as 



(4.1) 



This quantity can be expressed for paramagnetic case 
through the double occupancy and the number of d elec- 
trons as follows: S^_= Ud — 2D. 

Following Ref. [g^ we can distinguish two different 
regimes according to the values of the U Hubbard re- 
pulsion and the charge-transfer energy A. In the region 
U > A the insulating properties of the system are char- 
acterized by a charge-transfer gap. Instead, in the region 



U < A the system is in a Mott-Hubbard regime. Both 
regions are studied in the following. 

In Fig. and Fig. ^ we plot the squared local magnetic 
moment against the parameters A and U, respectively. 
5^ is an increasing function of of these parameters. Fig- 
ure |l| shows that the local magnetic moment takes 
the smallest value when A approaches zero; in such case 
p and d upper Hubbard levels coincide and due to strong 
mixing the double occupancy takes a large value. 
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FIG. 1. The squared local magnetic moment against the 
parameter A; the values of external parameters U, ut and T 
are given in the figure. The solid line is COM prediction with 
Pauli principle, the dotted one with decoupling. The squares 
are QMC data H. 
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FIG. 2. The squared local magnetic moment against the 
parameter U; the values of external parameters ut, A and T 
are given in the figure. The solid line is COM prediction with 
Pauli principle, the dotted one with decoupling. The squares 
are QMC data ||. 

Instead, when A becomes larger than U the system 
changes from charge-transfer insulator to Mott-Hubbard 
insulator and the local magnetic moment S'^ becomes 
independent on A. In this regime we are addressing a 
single-band 2D Hubbard model solution. 

An increasing behaviour of S"^ is also possible for 
charge-transfer insulator as it is shown in Fig. ^, where 
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A is taken equal to 2(7/3. In such case double occupancy 
of d electrons decreases with increasing distance between 
bands and the local magnetic moment saturates. 



La2-xSrxCuOi [Q, where it was stated that doping de- 
stroys antiferromagnetic spin correlations but does not 
destroy local magnetic moments. 
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FIG. 3. The squared local magnetic moment against tit 
for A = 1 and A = 4; the values of external parameters U 
and T are given in the figure. The solid line is COM pre- 
diction with Pauli principle, the dotted one with decoupling. 
The squares and the triangles are QMC data p9|. 



FIG. 5. Ud and Up against the total electron occupation 
TiT = rid + 2np for the values of external parameters [/ = 12, 
A = 2 and T = 1/8. The squares are Lanczos diagonalization 
results IBOl, the solid and dotted lines are COM calculation. 
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COM-d 
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FIG. 4. The tit-/! plot for the values of external parame- 
ters 17 = 6, A = 1 and T = 0.1. The squares are QMC data 
[ p9| , the solid line is COM prediction with Pauli principle, the 
dotted one is with decoupling. 

The dependence of the squared local magnetic moment 
on the total number of particles ut for A = 1 and A = 4 
and for the values of external parameters C/ = 6, T = 0.3 
is shown in Fig. |^. 

We also report results for the case when the parameters 
ds , fi X s are expressed by means of decoupling equations 
(^T9|)-(^^ (dotted hne). 

In all cases the agreement with QMC results is very 
good, specially in the case when we take into account 
the Pauli principle. 

From Fig. || we see that the local magnetic mo- 
ment only slightly changes with doping. Such fact 
was also observed in neutron scattering experiments on 



FIG. 6. Ud and Up against the total electron occupation 
riT ~ nd + 2np for the values of external parameters (7 = 6, 
A = 2 and T — 1/8. The squares and the triangles are 
QMC and Lanczos diagonalization results, respectively [ po[ ; 
the solid and dotted lines are COM predictions. 

In a previous paper the transition between p level 
and the lower Hubbard level ^ was not taken into account, 
and the the number of d electrons was calculated only 
approximately. As a result the agreement with QMC 
data was only qualitative. 

In Fig. ^ the dependence of ut on the chemical po- 
tential is drawn for [/ = 6, A = I and T = O.I; a good 
agreement between COM predictions and QMC results is 
found. As a comparison, it is also reported the case when 
the Pauli principle is not taken into account; in this case 
the agreement with QMC simulations is not so good. 

In Fig. 1^ and Fig. ^ we present also a very good agree- 
ment of our results with Lanczos diagonalization and 
QMC calculations given in Ref . pO] . The authors of this 
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work write the full number of holes as n = ncu + 2no, 
where the double number of oxygen holes is due to the 
two oxygen ions in the cell. We remind that in the p-d 
model there are two bonding and two nonbonding elec- 
trons per site. To make a comparison with results of Ref. 

we have to take into account the Fermi occupation 
number o f the se two unbending orbitals and to rewrite 
equation ( 0.1^ ) as 



riT = 4: - 2{Cii + C22 + C33) 
2 



exp((ep-/x)/T) + l- 



(4.2) 





■ QMC 

COM 

U=6 

T=1/8 











For the values of external parameters given in Fig. ^ 
and Fig. ^ the system is in the charge-transfer regime and 
in the vicinity of ut = 1 holes prefer to go to p orbitals 
and electrons to d orbitals. 

In Fig. I?] it is shown the dependence of the charge- 
transfer susceptibility at q = on charge-transfer energy 
A. This quantity |]30|| is given by 



lim xcT(q) 



_d_ 
dA 



{ud - 2np). 



(4.3) 



The agreement between COM predictions and QMC 
simulations is very good. 




FIG. 7. The charge-transfer susceptibility xct at q = 
versus A for the values of external parameters U — 6, ut ~ 5 
and T — 1/8; the squares are QMC results and the solid 
line is COM prediction. 

In Fig. Fig. || and Fig. |l^ we plot the quantities Ud, 
Si and D against the charge-transfer gap A for the values 
of external parameters U — 6, ut — ^ and T — 1/8. The 
agreement between COM prediction and QMC simula- 
tion is excellent. The dependence of the calculated 
quantities on external parameters has the same behaviour 
as in Fig. |^. We want to note once more that the holes 
introduced in Cu02 planes reside primarily on oxygen 
sites, as it is well-known experimentally. 



FIC. 8. The n^-A plot for the values of external parame- 
ters (7 = 6, riT ~ 5 and T = 1/8; the squares are QMC results 
[pO| and the solid line is COM prediction. 
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FIG. 9. The squared magnetic moment against A for 
the values of external parameters U — 6, tit = 5 and T = 1/8; 
the squares are QMC results ^] and the solid line is COM 
result. 




FIG. 10. The double occupancy D versus A for the val- 
ues of external parameters C/ = 6, = 5 and T = 1/8; 
the squares are QMC results js^ and the solid line is COM 
prediction. 
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In Fig. |Tl| and Fig. [T^ the band dispersion of Zhang- The density of states (DOS) for p and d electrons is 

Rice singlet for tit = 2.5 and ut = 2.75 is shown. The respectively given by the following expressions 
squares are QMC data from |31| and the solid line is the 

COM result. The data differ significantly only for the /" 2 v-^ (")/ 

values of dispersion energies far from the chemical po- ^p\^) = J ^ ~ £^n(k)) , (4.4) 

tential. Ob "^^ 
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FIG. 11. The band structure for the Zhang-Rice singlet for 
riT = 2.5. The squares are QMC data js^, the solid line is 
COM prediction with Pauli principle, the dotted one is with 
decoupling. 



^ -0.5 





■ 








■ QMC 

COM 

COM dec 

U=6 

n^=2.75 

1=0.1 

A=4 





(0,0) (71,0) (ii,ji) (0,0) 



A^<iM= / rf'kf][4^'(k) + 24S)(k) 

+a('3^)(k)]<5(c.-i?,„(k)), (4.5) 

where the spectral functions a"^^ (k) and the energy 
bands En (k) can be calculated as indicated in Sec. |l[ 

In Figs. I^, |l4| and |l| we present the DOS of d elec- 
trons (solid line) and p electrons (dotted lines) for C/ = 6, 
A = 2, T = 0.001 and for m = 2.70, 2.85 and 2.90, re- 
spectively. The solid vertical line at w = indicates the 
chemical potential. 
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FIG. 13. p - DOS (dotted line) and d - DOS (solid line) 
for riT = 2.70. With respect to the chemical potential 
Ep = -1.613. 



FIG. 12. The band structure for the Zhang-Rice singlet for 
riT ~ 2.75. The squares are QMC data j^, the solid line is 
COM prediction with Pauli principle, the dotted one is with 
decoupling. 

Summarizing, in this section we have presented a de- 
tailed comparison of COM results with the available data 
by numerical simulation. It is worth noticing that the 
formulation is fully self-consistent and no adjustable pa- 
rameter is used. 



B. Density of states 
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In this paragraph and in the next one we present an 
analysis of the band structures of the model and we give 
a description of the relevant physics near the Fermi level. 



FIG. 14. p - DOS (dotted line) and d - DOS (solid line) 
for UT = 2.85. With respect to the chemical potential 
e„ = -1.945. 
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FIG. 15. p " DOS (dotted line) and d ~ DOS (solid line) 
for riT ~ 2.90. With respect to the chemical potential 
Ep = -2.174. 

As it can be seen from these pictures there are four 
bands. The lower Hubbard band located at a; « —9 rep- 
resents the (, excitations; this band is filled mainly by d 
electrons and it has a small fraction of p electrons. The 
upper Hubbard band, located at w « 3, comes from i] 
operator excitations; it has mainly the weight of d elec- 
trons. The band coming from the Ep atomic level is lo- 
cated around « — 3. It is filled by p electrons, but 
due to strong mixing it also contains a large fraction of 
d electrons. Finally, the Zhang-Rice singlet band is situ- 
ated around the chemical potential; this band has almost 
equal fraction of p and d electrons. 

We want to stress the importance of studying the 
Zha ng-R ice excitation as independent field in the initial 
set (pTsI). As shown in Ref. |24|, if one does not initially 
consider the Ps excitation, the information about p and 
d elec trons coupling is lost after taking averages ( [2.20 ), 
(2.21) and one does not get the band dispersion situated 
at the chemical potential when the total number of par- 
ticles is ut ~ 3. 

The singlet coupling of p and d electrons was investi- 
gated by means of COM in the earlier work of Matsumoto 
et al. together with one of us . They obtained the sin- 
glet excitation band at the chemical potential and studied 
how the density of states depends on the variation of the 
correlators. 

In this work we considered the full p-d model, as de- 
scribed by Eq. (2.1), and we use d ad dit ional equations 
supphed by Pauh principle (Eqs. ( 3.15 ) -( [3.17| )) to calcu- 
late the values of these correlators. 



In Fig. |6| the band structure for the Zhang-Rice sin- 
glet is drawn for the values of external parameters given 
in the pictures; we see that at the momentum k = (tt, 0) 
such band has a saddle-point. This saddle-point leads 
to the van-Hove singularity in the density of states (see 
Fig. |l|, Fig;^ P and Fig. ^ and it was observed in ex- 
periments p-lQI. 

We note that for U = 6 and A = 2 the coincidence 
of the chemical potential and van-Hove singularity takes 
place at the value of total number of particles ut = 2.85, 
corresponding to the hole doping 6 = 0.15. This re- 
sults to an enhancement of the thermodynamical proper- 
ties such as specific heat and spin magnetic susceptibil- 
ity, as observed in the two-dimensional Hubbard model 
||,[|j3|. Upon increasing the particle number up to 
Ut = 2.9 the gap between the upper Hubbard band and 
the Zhang-Rice singlet band disappears. 

In Fig. |l^ and Fig. ^ we plot the Fermi Surface for 
the values of the external parameters U = 6, A = 2 and 
T = 0.001 and for total number of particles ut = 2.7 
and riT = 2.9 respectively. The values of the tempera- 
ture and of the total number of particles were chosen in 
order to compare the Fermi surface calculated by means 
of COM with the experimental one measured by Ino et 
al [ p4| and shown by the circles in the pictures. 

We see that at the total number of particles tit = 2.85, 
when the chemical potential crosses the van Hove sin- 
gularity, the Fermi surface changes its shape from the 
electron-like in the overdoped regime (n-r = 2.7) to the 
hole-like in the underdoped regime (nx = 2.9). This is in 
agreement with ARPES results |Q and Hall coefficient 
measurements p5[ ]. 

As a comparison, we also plot the Fermi Surface when 
we do not take into account the Pauli principle (dotted 
line) ; the agreement with experimental results is not good 
in this case. 



C. An analysis of the experimental data 
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We investigated the region of total number of electrons 
around ut = 3, which corresponds to the initial situation 

where only one electron is in the copper d orbital due to FIG. 16. The band structure for the Zhang-Rice singlet for 

the strong Hubbard repulsion and two electrons in the -"-t = 2.70, 2.85, 2.90 and for C/ = 6, A = 2 and T = 0.001. 
bounded p orbitals of oxygen. 
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FIG. 17. The Fermi Surface for f/ = 6, A = 2, nr = 2.7 
and T = 0.001; the sohd hne is COM prediction with Pauli 
principle, the dotted one with decoupling. The circles are 
experimental results pS. 
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FIG. 18. The Fermi Surface for ;7 = 6, A = 2, ut = 2.9 
and T = 0.001; the solid line is COM prediction with Pauli 
principle, the dotted one with decoupling. The circles are 
experimental results pS. 



V. CONCLUSIONS 

The p-d model was studied in order to describe ener- 
getic properties of charge excitations in the Cu02 planes 
of cuprates in the normal state. Excitations are described 
in terms of composite operators in the Cu02 cluster. The 
most important excitation appears to be the Zhang-Rice 
singlet. The oxygen hole doped into the Cu02 plane cou- 
ples with Cu electron forming a singlet state. The band 
of these states is located between the oxygen level and 
the copper upper Hubbard atomic level. By hole doping 
the Fermi level crosses this band and at definite value 
of doping it coincides with a van-Hove singularity in the 
density of states. This singularity is formed by the sad- 
dled shape of the band dispersion at the Brillouin zone 
points (tt, 0) and (0, tt) as it was also obtained in experi- 
ments and in QMC calculations. 

The large volume of the Fermi surface at low doping 
results from the small spectral weight of the field de- 
scribing p electrons dressed by spin-flips of d electrons. 
Such reduction of the spectral weight is possible due to 



its redistribution to the other bands strongly mixed with 
Zhang-Rice singlet band. 

By using a four-pole approximation in the framework 
of the Composite Operator Method various local quanti- 
ties have been calculated as functions of model and phys- 
ical parameters. The results of calculations have been 



presented in Sec. IVA and compared with the results of 
numerical analysis. 

We want to stress the relevance of Pauli principle 
in treating with strongly correlated electronic systems 
within our method. The comparison of COM results with 
numerical data, concerning local properties and band 
structures, and experimental measurements obtained by 
Ino et al looks to have a good agreement when we 
take into account the Pauli principle in computing the 
fermionic propagator. 

As last remark, we want to note that we can calculate 
the two-particle Green's functions in the one-loop ap- 
proximation by using the single-particle propagator. 
Calculations in this direction are now in progress. 
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APPENDIX: THE MATRICES I AND M AND 
THE SELF-CONSISTENCY PARAMETERS. 

In this appendix we give the explicit form of matri- 
ces / and m and we also report the expression of the 
self-consistency parameters. 

The / matrix is diagonal and its elements are: 



111 = 1, h 



1 - nd/2, I33 = nd/2, 



I44 = 3(nrf -2D - a{k)xs) + 

-9(— + — ). 
-122 -133 



(Al) 



(A2) 



The quantities Ud and D are the particle number and 
the double occupancy per site of d electrons, respectively: 

nd^{dHi)d{i)), (A3) 

D^{n^{t)n^{i)). (A4) 
The other parameters are defined by 

b^{p'y{{)r^Hz)), c=(p^(j)et(j)), (A5) 



as = (p^(i)ri(«)) , Xs^l K(0<(«)) • (A6) 

The m matrix is hcrmitician and its elements have the 
expressions: 
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mii=ep-/x, mi2 = 2i/227(k), 

mi3 = 2t/337(k), mi4 = 0, 

= (£d - m)-?^22 + 2t(c - 6), = 2(6 - c), 

m24 = 2t/7rp^, m34 = -TO24, 



"^33 = (Ed + - /^)^33 - 2i(6 - c), 
77144 — {£p — A^)^44 + '^tlksPs + ^p-^Jrp 



where 



-^^Ps = -^{rid-W - a(k)xs) + 2as - 



tp{b - c) 



2t 



Ik,^ = -^C^tlirp, ~ £pih2 - tp{c - b)), 



h,r) = --^C^tlTTp, + Sprihz +tp{b- c)), 



6(-5-/ + 4a(k)) -46, 



-3 I c— ^ + b——^ 

J22 -(ss 



The parameters fe,., /, and c?, are defined by 

/ = (p''(z)dt(z)p^(z)p'^'(z)V 



ds^-{nl{i)akp'{i)d^{i)). 



(A8) 

(A9) 
(AlO) 
(All) 

(A12) 

(A13) 
(A14) 

(A15) 
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